States of the Dirac equation in confining potentials. 
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We study the Dirac equation in confining potentials with pure vector coupling, proving the exis- 
tence of metastable states with longer and longer lifetimes as the non-relativistic limit is approached 
and eventually merging with continuity into the Schrodinger bound states. We believe that the ex- 
istence of these states could be relevant in high energy model construction and in understanding 
possible resonant scattering effects in systems like Graphene. We present numerical results for the 
linear and the harmonic cases and we show that the the density of the states of the continuous spec- 
trum is well described by a sum of Breit-Wigner lines. The width of the line with lowest positive 
energy, as expected, reproduces very well the Schwinger pair production rate for a linear potential: 
we thus suggest a different way of obtaining informations on the pair production in unbounded, non 
uniform electric fields, where very little is known. 
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The absence of bound states for the Dirac equation in 
confining potentials poses a delicate question of physical 
interpretation. It was in fact shown in [l[ that the asymp- 
totically oscillating behavior of the solutions of the Dirac 
equation with vector coupling and a potential given by 
a positive power of the modulus of the position variable 
implies a purely continuous spectrum. The result has 
been subsequently confirmed and generalized by prov- 
ing that even any self-adjoint extension of the boundary 
value problem has only a purely continuous spectrum, [2| . 
This situation is contrary to the physical intuition and 
it makes hard to justify a perturbative approach to the 
relativistic corrections, since no definite perturbed eigen- 
values exist. Most of later investigations have therefore 
tried different ways to introduce confining potentials into 
the Dirac equation, e.g. by scalar coupling or by pro- 
jection onto the large component (see [3| for a review): 
doing like that, however, physically relevant systems as 
charged particles in strong electric fields do not find an 
appropriate description. A possible way out of this dif- 
ficulty was given in the mathematical paper 0], where 
a so called 'weak quantization' was introduced to treat 
the (l+l)-dimensional Dirac equation with a linear po- 
tential. The analysis is developed in the complex plane 
of the energy; when the Schrodinger limit is approached 
it is shown that the real part of the energy converges 
to the non-relativistic spectrum and the imaginary part 
becomes exponentially vanishing. For a linear potential 
analytical solutions are available in terms of special func- 
tions and in 0] the spectral quantities were estimated by 
perturbative expansions, having thus a very limited range 
of validity, out of which no general picture can emerge. 
Moreover a coherent physical interpretation was clearly 
outside the mathematical purpose of the author. To our 
knowledge no further development along these lines is 
found in literature since. 



From a physical point of view, the presence of an un- 
bounded increasing potential brings to bear upon the 
problem arguments similar to those of the Klein paradox, 
widely studied both in first and in second quantization, 
[H, El 0] • Recently a field-theoretical interpretation using 
numerical methods based on spatial and temporal reso- 
lution was given in where it is found that the pair 
production by the potential is suppressed when the spa- 
tial density of the incoming electron overlaps with the 
potential region and that the transmitted portion of the 
wave packet, in a single particle description, corresponds 
to the amount by which the electron reduces positron's 
spatial density. Although in this letter we strictly re- 
main in a first quantized framework and we are dealing 
with the stationary problem defined by the Dirac equa- 
tion rather than with a scattering picture, still we shall 
see that the pair production rate is recovered in a natural 
way. Our treatment follows the classical methods of the 
spectral analysis, [9l.ll0l|. and, in particular, it deals with 
an accurate evaluation of the density of the states of the 
continuous spectrum. It can be compared with the phe - 
nomenological approach in terms of Gamow vectors, [llj . 
mainly used to describe resonances in composite systems 
of solid state and nuclear physics. The investigation is 
necessarily numerical: thus it does not suffer the limi- 
tations of the perturbative expansion and it can easily 
be extended to more general potentials for which ana- 
lytical solutions do not exist. We will present in detail 
the results of the linear and quadratic potentials. The 
general situation can be summarized as follows. In non- 
relativistic quantum mechanics the spectrum is discrete, 
the eigenvalues correspond to the real zeroes and poles Ai 
of a spectral function introduced by Weyl (traditionally 
denoted by ra(A), [§]) and the density of the states reduces 
to a sum of <5-functions, one for each bound state. In the 
Dirac equation the A^ move off the real axis, the spectrum 
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becomes purely continuous and the density of the states, 
p(A) = — 3?(m(A))+3?(l/m(A)) for real A, appears now as 
a sum of Breit-Wigner (BW) lines whose central values, 
determined by 5ft(A;), are closer and closer to the non- 
relativistic eigenvalues and whose widths, determined by 
S(Ai), are more and more narrow for decreasing values 
of the ratio of the interaction to the mass energy. 

The physical interpretation suggested by these facts is 
that the broadening of the 5-lines is due to transitions be- 
tween positive and negative energy sectors induced by the 
supercritical field: thus, although in the relativistic con- 
text the nature of the spectrum is completely changed, 
still narrow BW lines signify the presence of metastable 
states, giving rise to resonances in the scattering cross 
section around the line energies. On the one hand, there- 
fore, the continuity to the non-relativistic states is re- 
covered. On the other, contrary to what occurs for the 
Schrodinger ground state, in the relativistic regime also 
the lowest positive energy state decays. The second quan- 
tized counterpart of this fact is that the Fock vacuum will 
not remain such forever, but, according to the usual the- 
ory of the effective action, [o, EH , it will decay with the 
exponential law |(0(t)|0(i + T))| 2 = exp( 
w-f is the pair production rate for unit volume and unit 
time. We thus expect similar behaviors of the line width 
of the lowest positive energy state and of the pair pro- 
duction vs. the interaction strength. This circumstance 
appears very well verified for the pair production rate in a 
constant electric field, as obtained by Schwinger 1^, 16| , 
so that we are led to assume that the width of the first 
resonance can provide a quantum mechanical way of esti- 
mating the pair production for general situations where 
little is known from QED, as in the case of unlimited 
growing potentials (see H, 14 1 for recent developments) . 
We then present new data for a quadratic potential, cor- 
responding to an electron in a uniformly growing electric 
field, finding a pair production behavior very similar to 
the Schwinger's one. We finally believe that our results 
can be relevant not only in model building of quark sys- 
tems, 17, [III], but also in investigations of Klein para- 
dox in strong crystalline fields [l9[ as well as in the very 
recent and expanding subject of the Graphene physics, 
where the influence of impurities is described by the Dirac 
equation with vector coupling, [2(|: the metastable states 
may result essential for understanding possible effects of 
resonant scattering. 

Consider the (l+l)-dim Dirac equation in a unit sys- 
tem with h = 1 , 

iP'{x) - [ (l/c)(E - U{x)) io v +mca x ] 1>(x) = (1) 

where ip(x) = t (il>i(x), ip2 {x)) and cr, are the Pauli ma- 
trices. For the family of even potentials U(x) = a\x\ n 
the equation ([1]) can be studied in [0, oo), having in- 
finity as the unique singularity in the limit point case 
Therefore, from the Weyl general theory of singular 



boundary value problems, f° r ^(E) > there ex- 

ists only one normalizable solution ip(x) of the equation, 
up to a constant factor. If {tp^ (x, E)} i= i^ is a funda- 
mental system of solutions of (fTJ) with %l>y{Q,E) = 5j, 
the Weyl function m(E) is defined by an expansion 
if(x,E) =tpW(x ) E)+rn(E)ipW(x,E). Using the initial 
conditions, m(E) can be written as 

m(E) = M0,E)/M0,E) (2) 

Finally the density of the states for real Eq reads, (Toj . 

p{Eo) = Hm (-9(m(J^) + ii/)) + 9(l/m(Eb + w))) (3) 

For calculation reasons, we find it convenient to define 

<t>{x) = t (0i(x),0 2 (a;)) =2-ii(a v + <T z )1>(x) (4) 
obtaining for 4>{x) the following equation: 

<P'(x)- [(i/c)(E - U(x)) a z - mca x ] <j)(x) =0 (5) 

We denote by {<f>^(x,E)}i=i,2 the fundamental solu- 
^Tw^)" where tions of ([5]) with initial conditions (f>^ (0, E) = 5j . Hence 



if 4>{x 1 E) is the normalizable solution of ([5]) correspond- 
ing to ip(x,E) and we expand <f)(x,E) = (^'(x^E) + 
k{E)4>( 2 \x, E), we then have a finite limit 

k(E) = - lim <fP{x,E)/(j)f\x,E), i = 1,2. (6) 

x — >oo 

The relation between m(E) and k(E) is given by 

m(E) = i (k(E) + i) (k{E) - i)' 1 (7) 

Finally, introducing the "non-relativistic energy" E Nn by 
E = £ NR +mc 2 , the elimination offeix) yields the second 
order equation for <j>\ (x) 

4>'{{x) + [2m(E mt - U(x))+c- 2 R{x)\ ^{x) = (8) 

where R(x) = icU'(x) + (E NR - U(x)) 2 . The non- 
relativistic limit for c — > oo is then evident. 

Let us now consider the specific cases of the potentials 
U(x) — £\x\ and U(x) = (l/2)muV. For the linear 
potential we introduce 

y = (2m£)ix, X=(^E m , ^ = ^(^) § (9) 

and for the quadratic potential we let 

V = (mw)3 x A = (2/w) E NR , n = W /(4mc 2 ), (10) 

so that in both cases the non-relativistic limit is obtained 
for fi — ► 0. Equations |T|), (O and ((8|) , then, specify to 

Qi[A n (y)ia v + (2Q)- 1 a x ]ij(y) = (11) 

4>'{y) - i & [A n {y)a z + i(2Q)- x a x ] 4>{y) = (12) 

4>'({y) + [nA 2 n (y) + l niny n - 1 - Mv) = (13) 
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FIG. 1: The density of the states for the linear potential with 
= 0.3, 1. The scale for Q, = 0.3 must be multiplied by 10 2 . 
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FIG. 2: p(A, fi) for the quadratic potential. 



with n = 1, 2 and A n (y) = A + 1/(20) - y™. When n = 1 
the normalizable solutions of (|13p with complex spectral 
parameter are known, [4|, and they are all proportional 
to the cylinder function Di T (—z), with r = (2J1 1 / 2 ) -3 
and z = (— 4Sl) 1 ' 4 Ai(y). Carrying out the calculations 
previously described, it is straightforward to arrive to the 
expression for the density of the states for equation pip . 
The general properties of p(A, f2), as in ([3]), can be appre- 
ciated by looking at the complete numerical results. First 
we see the convergence to the Schrodinger levels when 
— > 0. For instance, the peaks of the first resonances for 
fl = 0.01 are located at 1.0197, 2.3274, 3.2284, 4.0555, 
4.7756: these should be compared with the first zeroes of 
the derivative of the Airy function and of the function it- 
self (even and odd solutions respectively), 1.0190, 2.3384, 
3.2482, 4.0884, 4.8201. Secondly, in comparison with the 
non-relativistic case, the spacing of corresponding reso- 
nances decreases for increasing values of the energy and 
of the interaction strength. In FigQ]we present the plot of 
p(A, fi) for n = 1, f2 = 0.3 and 1. The fit of the density of 

the states by a sum Yh=i a A w7i((^ 2 ''~^oif '+^lf1i)~ 1 of 
BW curves with the appropriate parameters Xoi, 7j and 
coefficients Cj, gives a perfect superposition with p{\ 1 Q). 
In the left part of FigH] we compare the width 7 of the 
first resonance with the pair production per unit length 
and time in a uniform electric field, [15j, that in the vari- 
ables (0 reads 



In 



1 — exp 



(-7r/(4f2 3 / 2 ))] (14) 



As stated above, the excellent agreement, without any 
free parameter to be adjusted, proves that pair produc- 
tion and line broadening are two different descriptions of 
a same physical situation. The minor differences, mainly 
for increasing f2, can be partly assigned to the description 
of the spectrum in terms of BW lines and partly to the 
fact that the Schwinger pair production is an effective 
one-loop calculation, possibly under-estimating the ac- 
tual production rate [13j : one could make these small dif- 
ferences vanishing not by exponential but only by power 



law corrections in f2. 

Similar considerations also apply to the quadratic po- 
tential (n = 2 in (|1HI13|) ), although the computational 
technique is now different. Actually, solutions of (113p 
exist in terms of triconfluent Heun functions Ht, [2l|, 
(j>i(y) = Ae f H T (p, -3,<?, z)+Be^ i H T (p,3,q,-z) , where 
/ = (A + (2Q)- 1 - y 2 /3) y, p = (3 (4fi)- 2 )l , q = 

(12 ft)i (A + (2n)- 1 ) and z = -i(2Qs/3)7 y. Unfortu- 
nately no sufficient informations on the asymptotic be- 
havior of Ht are available, to our knowledge, to deter- 
mine the normalizable solutions for complex A. Hence we 
use a completely numerical scheme that extends to any 
potential U(x) = a|x|™, for which analytical solutions do 
not exist when n > 3. 

The calculations are straightforward and follow step by 
step the theory we have previously summarized. First we 
find, by numerical integration, a fundamental system of 
solutions of (|12p . from which we determine k(X) accord- 
ing to ([S]). Some care must be used in taking the limit, 
that is approached not in a monotonic but in an oscil- 
lating way, as it is evident by looking at the asymptotic 
leading terms of (fl3|) : the convergence is increased by 
constructing the sequence of the average points of pairs 
of nearby maxima and minima, whose limit is searched 
with sufficiently high absolute and relative precision. We 
then find the Weyl function m(A) from ([7]) and eventually, 
from ([3]), we deduce the density of the states, see Figl2] 
The maxima of the first four BW lines are displayed in 
FigEl Starting from the odd integers, that correspond 
to the non-relativistic values, we see that their spacing 
decreases both for increasing and A: the same effect 
has been observed for the relativistic Landau levels, [13] , 
and it loosely seems to repropose, in relativistic quan- 
tum mechanics, the usual relationship between circular 
and harmonic motion. We can also remark that, as for 
the linear potential, the lowest resonance has a central 
value that always remains near the non-relativistic value. 
The data of the higher resonances, instead, are well fit- 
ted by decreasing exponentials in f2, they approach to 
each other and their unequal spacing should be taken 
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FIG. 3: The first BW maxima vs. Q for the quadratic poten- 
tial. Non-relativistic bound states correspond to odd integers. 
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FIG. 4: Left plot: the width of the first resonance (diamonds) 
compared with the pair production curve ur (fi) (solid line). 
Right plot: the width of the first resonance (circles) for a 
quadratic potential. The solid line, giving a very good fit, is 
w(Q) = -7r _1 ln[l -exp(-7r/(4 2 fi 3/2 ))], analogous to w f (Q.). 



into account in the construction of quarkonium models. 
The right part of Fig Q] reproduces the width of the first 
BW resonance vs. ft. According to what we said above, 
the plot can give an estimate of the pair production for 
this case, not treated by QED. Remark that the data are 
well approximated by a curve w(il) very similar to the 
production in a constant field. 

In conclusion, by using the possibilities offered by 
quantum mechanics, we have proved the continuity from 
the non-relativistic bound states to the states of the Dirac 
equation for an entire class of confining potentials, stress- 
ing the fundamental role of the density of the states and 
clarifying the apparent physical contradictions of the ab- 
sence of bound states for such potentials. Finally, al- 
though it is well known that the physical interpretation 
of the Dirac equation as a relativistic one particle wave 
equation presents many subtleties, we still believe that 
informations about its solutions and their properties de- 
serve a detailed knowledge. 
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